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Abstract 

We introduce a system of phenomenological equations for Bose-Einstein condensates of magnons 
in the one-dimensional setting. The nonlinearly coupled equations, written for amplitudes of 
the right-and left-traveling waves, combine basic features of the Gross-Pitaevskii and complex 
Ginzburg-Landau models. They include localized source terms, to represent the microwave 
magnon-pumping field. With the source represented by the <5-functions, we find analytical so- 
lutions for symmetric localized states of the magnon condensates. We also predict the existence 
of asymmetric states with unequal amplitudes of the two components. Numerical simulations 
demonstrate that all analytically found solutions are stable. With the 5-function terms replaced 
by broader sources, the simulations reveal a transition from the single-peak stationary symmetric 
states to multi-peak ones, generated by the modulational instability of extended nonlinear-wave 
patterns. In the simulations, symmetric initial conditions always converge to symmetric stationary 
patterns. On the other hand, asymmetric inputs may generate nonstationary asymmetric local- 
ized solutions, in the form of traveling or standing waves. Comparison with experimental results 
demonstrates that the phenomenological equations provide for a reasonably good model for the 
description of the spatiotemporal dynamics of magnon condensates. 

PACS numbers: 75.40. Gb; 03.75.Nt; 74.20.De 
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I. INTRODUCTION 



The experimental observation of the Bose-Einstein condensation (BEC) in dilute gases 
of alkali atoms [l| has been a milestone in the development of atomic and condensed-matter 
physics, demonstrating the reality of this state of matter, and providing a unique testbed 
for studying numerous macroscopic quantum effects [2j. As is well known, the transition 
to the condensation in atomic gases occurs at extremely low temperatures, ~ 10~ 7 K, the 
number of atoms in the condensate usually being quite small, < 10 4 . 

These achievements were followed by the creation of BEC in gases of bosonic quasi 



particles representing fundamental excitations in solid-state media, including excitons 



polaritons 



a 



a, 



4|, triplons [5| and magnons (spin waves in ferromagnets) [6|. Similar to the BEC 
of atomic gases, the condensation of quasi-particles takes place when their density exceeds 
a certain critical value, which increases with the temperature of the system. However, 
the advantage of studying BEC in solids is that the density may be increased by external 
pumping fields: laser or microwave radiation, in the case of excitons and polaritons, or 
magnons, respectively. Therefore, the BEC transition in solid-state media can be reached 
at much higher temperatures than in the atomic gases. 

Unique properties are featured by the BEC of magnons in tangentially magnetized 
yttrium-iron-garnet (YIG) films. In particular, the condensation of magnons occurs at room 
temperatures Moreover, while, in the first experiments, monochromatic pumping 

fields were used, it has been recently reported that the BEC of magnons can be achieved by 
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Despite the achievements in the experimental studies of BEC of magnons, development 

R ^, 111 Thermodynamics of the quasi- 
equilibrium magnon gas and the stability of the condensate were considered in Refs. js] 



means of an incoherent microwave pumping 
Despite the achievements in the experime 
of the corresponding theory is just starting 
lib] 

n 

and [9] respectively. The formation of the coherent magnon state (the BEC proper), accom- 
panied by the emergence of macroscopic dynamic magnetization, is described in Ref. jiol ]. 
The theory is capable of explaining, in a quantitative form, the observed intensities of the 
Brillouin light scattering and microwave signals generated by the condensate. The presented 
work addresses the dynamics of the magnon condensate, first experimental studies of which 
were recently reported in Ref. [l^ . 

The magnon BEC is fundamentally different from its counterparts in all other systems. 
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The condensation of magnons occurs at two separate points in the space of the magnon fre- 
quency and wavenumber, (uo, ±k ), hence, for the theoretical description of the condensate 
one needs to introduce wave functions with two components corresponding to these points. 
One should take into account interactions between these two components, which may result 
in effects which are not observed in BEC of other (quasi-)particles. 

In accordance with the above discussion, different types of models are necessary for the 
description of BEC in atomic gases and solid-state media have to be used. In the former 
case, a widely adopted theoretical approach is based on the Gross-Pitaevskii equation (GPE) 
for the single-atom wave function 2[]. In the mean- field approximation, the GPE takes into 
account collisions between atoms through the cubic term. In fact, the GPE provides for a 
very accurate model for the practically relevant case of dilute ultra-cold gases. An inherent 
property of the GPE is that the norm of the wave function is a dynamical invariant, which 
complies with the obvious condition of the conservation of the number of atoms in the 
gas. Equations that may serve as a mean-field model of the condensate of quasi-particles 
should be different, as the quasi-particles (in particular, magnons) may be created by the 
pump field and lost through linear and nonlinear dissipation. Therefore, the respective 
equation is expected to be of the complex-Ginzburg-Landau (CGL) type [13], combining 
conservative and dissipative terms. However, a difference from the usual form of the CGL 
equation is that, in the present case, it must include a source term representing the localized 
pump (in this respect, the equation may be somewhat similar to those describing nonlinear 
optical cavities pumped by external laser beams [14]). Nevertheless, this difference is not 
absolute, as ultra-cold atoms have an effective finite lifetime, due to heating effects and 
three-body interactions. As discussed in Refs. 1^, the ratio between the lifetime of the 
(quasi-)particles and their thermalization time is decisive for this issue. For magnons this 
factor takes values in the range of 5 — 10, whereas for ultra-cold alkali atoms it may reach 
two orders of magnitude. 

The objective of the present work is to put forward a semi-phenomenological system of 
coupled equations for the wave functions of the magnon condensate generated by the local 
microwave source. Comparing predictions of the model with experimental observations of the 
dynamics reported in 12|, we conclude that the equations correctly reproduce main features 
of the spatio-temporal dynamics (pattern formation in the condensate). The equations are 
introduced in Section II, and some exact analytical solutions for localized states, predicted 
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by the model, are reported in Section III. This is followed by the presentation of numerical 
solutions in Section IV and, eventually by the comparison with the experiment in Section 
V. The paper is concluded by Section VI. 



II. THE MODEL 

As mentioned above, the onset of the condensation in the magnon gas reveals itself 
in the emergence of the macroscopic dynamical magnetization in the sample. It should be 
emphasized that only condensed magnons, with the above-mentioned values of the frequency 
and wavenumbers, (uio,±ko), contribute to this magnetization, whereas all other magnons 
merely decrease the absolute value of the static magnetization. Following the envelope 
(slowly- varying-amplitude) approximation, widely adopted for the analysis of nonlinear spin 



waves in ferromagnetic films [la, ll7|, [18] , as well for the description of nonlinear waves in 



many other settings, we introduce the order parameter (the full wave function of magnons) 

as 

V(z,t) = ^ + {z,t)e lkoZ + V„(z,t)e~ ikoZ , (1) 

where absolute values of amplitudes represent dynamical macroscopic magnetizations 
created by the corresponding condensate components, m±, normalized to the absolute value 
of the static magnetization (Mo): \^±\ 2 = \m±\ 2 / (2Mq), while phases of ^± account for 
the coherence of the two components of the condensate. In Eq. (JTJ) , it is assumed that the 
magnetic field is aligned with axis z, the microwave antenna is oriented perpendicular to it, 
and ±fco are large wavenumbers (ko — 3 x 10 4 cm -1 in YIG films) at which the transition 
into the BEC state takes place. The remaining dependence on z in ty±(z,t) is assumed 
to be slow in comparison with the rapidly oscillating carrier waves, exp (±ikoz), therefore 
the interaction between and is incoherent [see Eqs. ([2]) and ([3]) below]. The slow 



dependence arises due 



In Refs 



. Q,[l7 



ocal variations of the pumping which excites magnons. 



18( | it was shown that, in the case of single-component magnon wave 
function, the evolution of the envelope is described by the nonlinear Schrodinger equation. In 
a similar way, for two components \]/± one arrives at a system of phenomenological equations 
of the GPE/CGL type, which include the above-mentioned source terms, written below as 
—fS(z) to allow analytical treatment of the problem (numerical results will be reported 



4 



below for a realistic broad shape of the source). In the scaled form, the equations are: 

i (*+) t + Sfio ■ ^+ + - {^+) zz + ir)*+ 
+ (|^ + | 2 + a 1 |^_| 2 ) m + + it (|^ + | 2 + o 2 |^_| 2 ) ^+ = -/ 5{z). (2) 



1 

2 

(|^_| 2 + (7! |^+| 2 ) ^_ + *r (|^_| 2 + a 2 |^+| 2 ) = -/ 5(z). (3) 



Here, <5/io is a possible shift of the chemical potential of magnons in the condensate with 
respect to that of uncondensed magnons. We also introduce phenomenological damping, 
represented by linear and cubic terms in Eqs. (j2D and (jHD, with coefficients r\ and r, 



respectively. The same nonlinear-damping term was introduced in Ref. [19|. Further, <Ji, o 2 
are the XPM/SPM (cross/self-phase-modulation) ratios for the conservative and dissipative 
nonlinear parts of the equations. The sign in front of the SPM and XPM nonlinearities 
in the conservative part corresponds to the case of attractive magnon-magnon interaction, 

nn 

which was realized in tangentially magnetized YIG films in experimental works [9|, H8( • 

It is relevant to mention that equations similar to Eqs. Eqs. (TSJ) and ([3]) were derived 
by Lvov for the description of spatially inhomogeneous parametric excitation of magnons 



20| | (in that case, two groups of magnons with wavenumbers ±k were also excited). In 
fact, the applicability of amplitude equations of this general (Ginzburg-Landau) type to the 
description of condensates is a universal fact (at least, at the phenomenological level), even 
in the case of strong interactions between quasi-particles 21| . 

The notation adopted in the equations implies that t is normalized to the characteristic 
time of the nonlinear interactions of magnons, r ~ 3.6 ns, which, as said above, is much 



shorter than the lifetime of a magnon (T m ~ 250 ns in YIG [221]). hence rj = r / (2T m ) is 
considered as a small parameter. Spatial coordinate z is normalized to characteristic length 
zq, which is defined by the ratio of the nonlinear interaction strength to the dispersion 



coefficient (zq ~ 0.5 /im. in YIG |18|). Unlike other coefficients, which may be identified 
from empirical data, 5{Iq should be treated as a phenomenological parameter, that may be 
found from adjusting theoretical predictions to experimental observations. 

The strength of the source, /, may be set to be real and positive, by definition. As 
mentioned above, function 5(z) in Eqs. ([2j), (j3J) is the Dirac's 5-function if we aim to find 
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analytical solutions (see below), or a regularized counterpart of the 5-function for numerical 
solutions. In the latter case, it is taken as 

A being the regularization parameter. 

Note that Eqs. ([2j), (j3J) are somewhat similar to the coupled CGL equations for counter- 
propagating waves in the binary- fluid thermal convection [23[ . However, that model includes 
finite group velocities, but not the source terms. Equations (TSJ), (jHJ) do not contain group- 
velocity terms [which might be, generally speaking, ±ic (\&±) 2 ], because the group velocity 
of magnons vanishes at the condensation points, k = ±/co- 

In the case of the 5-function, symmetric stationary solutions to Eq. ([1]) may be sought for 
as ^+(z, t) = *&-(z) = ^(z), where the single complex function, *&(z), satisfies an ordinary 
differential equation, 

(5fM> + irj)^f+ ]^ + [(1 + at) +tr(l + a 2 )\ |tf | 2 tf = 0, (5) 

at z / 0, supplemented by the boundary condition (b.c.) at z — 0, which is generated by 
the integration of Eqs. (j2J) and ([3]) in an infinitesimal vicinity of z = 0: 

*'|*=+o - *Wo = -2/- (6) 

The meaning of this b.c. is that function ^(z) must be continuous at z = 0, but its 
derivative makes a jump at the position of the source. This approach to searching for 
stationary solutions is similar to that which was recently developed for a CGL equation 
(with an intrinsic gain, but without source terms) driven by the delta-like gain [241 ] . 

Concluding this Section, we emphasize that the source term introduced in the above 
equations is not direct proportional to the microwave field used in the experiment for the 
parametric pumping of magnons. In fact, the parametric pumping creates primary magnons 
away from the condensation points in the phase space. The BEC condensate then forms as 
a result of the thermalization of the primary magnons. Thus, the proposed model does not 
aim to explicitly describe the spatial spreading of the magnons during the thermalization 
process. 
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III. ANALYTICAL RESULTS 



A. Symmetric solutions 



Following the approach of Ref. [24j, where particular exact solutions were found in the 
framework of the above-mentioned CGL equation with the delta-like linear gain, we look for 
a localized solution to Eq. ([5]) in the following form: 

Ae ix 



(7) 



[smh(K(\z\+m 1 ~ iU, 

with real constants A, x, k, v and £. In particular, £ must be positive (otherwise, the ex- 
pression gives rise to a singularity at \z\ = — £). Actually, this ansatz was suggested by a 
formal singular solution to the cubic CGL equation with constant coefficients, that^Jn turn, 
is a counterpart to the so-called Pereira-Stenflo soliton, with sinh replaced by cosh 25|. The 
latter is a non-singular but unstable exact solution, which is valid in the case of rj < (a 
uniform linear gain, instead of the loss). Solutions of the Pereira-Stenflo type may be made 
stable as exact solutions to coupled CGL equations one of which is cubic, and the other one 
linear 261 ] . 

The substitution of ansatz (|7j) in Eqs. (|SJ) leads to the following equations: 



k 2 (2 



3iu) + 2 [(1 + m) +ir (1 + a 2 )] A 2 = 0, 



k 2 (1 - iv) 2 + 2 (fl + IT)) 



0. 



(8) 
(9) 



Further, inserting ansatz (J7]) into b.c. ([6]), we arrive at an additional relation, 

cos1i(k£) fe~ lx 



[sinh 



- zV) 



(10) 



Equations (jHJ) and (jHJ) can be solved directly [actually, it is the single physically relevant 
solution corresponding to ansatz ([7j)]: 



v 



, 9(l + ai) 2 | 3(1+^) 



^ = 5 

i/ ~ 2 



4(l + a|)r 2 2(l + <7 2 )r ' 



'2 + 



(i + ^ir 

1 + (T2 )V (l + a 2 )r 
2r (1 + a 2 ) ' 



(11) 
(12) 

(13) 



6n> = \iy 2 - 1) k\ (14) 

It is worthy to emphasize the meaning of Eq. ()14p : the solution given by ansatz (J7J) 
exists solely for the particular value of <5/i given by this expression. For other values of 
Sfjio, localized solutions are also expected to exist (see numerical results reported below), 
but they cannot be found in the exact form corresponding to the present ansatz. 

Once constants u, k and A are known, as per Eqs. (lllj) - ( fl3l) . the equation corresponding 
to b.c. (flOj) can be solved analytically too, and it also produces a single solution, which 
determines the remaining arbitrary constants, £ and x- 



cosh«) = ^/1 + — + — , (15) 

X = arctan {y) — z/ln (sinh (/«£)) . (16) 



A noteworthy feature of solution (TTojl is the absence of a threshold: any pump strength /, 
even a very small one, supports the solution. In the limit of / — > 0, Eq. ( fl~5l) yields £ — > oo, 
which naturally implies the disappearance of the BEC state when the microwave source is 
switched off. 



B. Asymmetric solutions 

A challenging issue for the theory is a possibility of the existence of asymmetric solutions, 
with ^+(z) 7^ ty-(z). In that case, the stationary version of Eqs. (j2J) and (jHJ) takes the 
following form: 

+ (|^+| 2 + a x |^_| 2 ) ^ + + it (|^ + | 2 + a 2 |^_| 2 ) ^+ = -/ 8{z). (17) 

+ (|^_| 2 + cti |^+| 2 ) ^_ + ir (|^-| 2 + a 2 |^+| 2 ) ^_ = -/ 5{z). (18) 

Here, we assume that S(z) is not an infinitely narrow delta-function, but rather a regularized 
one - e.g., as given by Eq. We will demonstrate below that in this case an asymmetric 
solution appears spatially nonuniformly, i.e., in regions of z where 5(z) exceeds a certain 
threshold. 
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The possibility of the existence of asymmetric solutions can be investigated analytically 
under the assumption that the linear and nonlinear dissipative coefficients, rj and r, are large 
parameters. Strictly speaking, this assumption contradicts the above-mentioned condition 
for the applicability of the equations to the description of the physically relevant situation, 
T] <C 1. Nevertheless, the analysis makes sense, in view of the fundamental significance the 
possibility of the symmetry breaking in models of the present type. 

Assuming that r\ and r are large, In the lowest approximation we neglect all conservative 
terms and replace Eqs. (ITT)) and (ITS)) by the following ones, for \P±(z) = i$±(z), where 
functions $±(z) are real: 

r ? $ + + r($ 2 + + a 2 $ 2 )$ + = /^). (19) 

v <5>_+t(<S> 2 _ + ct2$1)$- = f 5(z). (20) 

If an asymmetric solution exists, taking the difference of Eqs. (IT9l and (1201) and dividing it 
by A$ = $ + — $_ 7^ yield the following equation: 

(02 - 1) - ($ 2 + + $ 2 _) = r//r. (21) 

At the bifurcation point, where the asymmetric solution splits of from the symmetric one, 
$+ = $_ = $symm, Eqs. ( (Tl9l) . (120)) and (J2TJ) must be satisfied simultaneously, i.e., 

r/$symm + T (1 + 2 ) ^ symm = f 6(z), 

(22) 

(02 - 3) $ s 2 ymm = T]/T. 

The second equation in system (l22l implies that the bifurcation is only possible if 02 is 
large enough, viz., a 2 > 3. Then, the substitution of the second equation into the first one 
determines the bifurcation point, 

Further, equating <S>^±, as given by Eq. (123)) and by the second equation in system (122)) . 
yields the value of 5(z) at the bifurcation point: 



The meaning of Eq. (I24I) is that it gives a threshold condition for the symmetry-breaking 
bifurcation. Namely, the bifurcation occurs if the maximum value of S(z)f, which is its value 
at the central point, z = [for instance, fS(0) = f / (a/7tA) in the case of 5{z) given by Eq. 
(QJ] exceeds the threshold value: 



Mo) 



> 



Mo) 



(<7 2 -lf 47J 3 



(25) 



'- threshold (o~2 — 3) T 

In the present approximation, which corresponds to Eqs. f|T9|) . (120]) . that do not contain 
derivatives and therefore admit solutions with discontinuities of the first derivatives, the 
solution is asymmetric in the region of \z\ < Zut, where zui is to be found from Eq. ( 1241) . 
and the solution continues as the symmetric one to \z\ > zuf- This structure of the solutions 
resembles phase-separated states in two-component atomic BEC, found from the respective 



system of coupled GPEs in the Thomas-Fermi approximation 



27|]. 



IV. NUMERICAL RESULTS 
A. Stable symmetric solutions 

The numerical analysis was performed by means of direct simulations of Eqs. (T5]), (J3j), us- 
ing the Crank-Nicolson scheme. Starting with symmetric initial conditions, the simulations 
of the equations, where the regularized delta-function was approximated by Eq. (jlj) with 
appropriate values of A (in any case, it must be much smaller than the width of the estab- 
lished state), would always drive the solution to a stable symmetric localized state. Varying 
symmetric initial conditions for fixed parameters of the equations, we always observed the 
convergence of the numerical solution to the same stationary state, which is clearly an at- 
tractor of the model - at least, in the class of symmetric initial conditions (see a discussion 
of asymmetric configurations below). 

In all cases when parameter SfiQ was chosen as per Eq. (TH|) . so as to support the existence 
of the analytical symmetric solution found above, the established state produced by the 
simulations with A < 1.5 virtually exactly coincides with the analytically predicted form 
(thus confirming the stability analytical solutions). Typical examples are shown in Fig. (TJ for 
two different values of parameter / which accounts for the strength of the pumping. Other 
constants used for these simulations correspond, after all rescalings, to physical parameters 
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FIG. 1: Analytically predicted and numerically generated profiles of the localized states supported 
by the <5-function source, at parameter values = 0, ry = 0.007, r = 0.22, and u\ = oi = 2, A = 1 
for two different values of the source's strength, / (the small value of rj used here agrees with the 
condition of the applicability of the phenomenological model, see the text). 

relevant to the experiment (in particular, o\ = 2 is the value of the XPM/SPM ratio 
appropriate to the nonlinear wave interactions mediated by the cubic terms |28j). 

Stationary symmetric solutions of another type were produced by the simulations in the 
case of a broad source, represented by Eq. (00) with much larger values of A: the increase of 
the source's strength, /, and the subsequent growth of the amplitude of the stationary state 
leads to a transition from the simple localized shape with a single maximum to multi-peak 
patterns, as illustrated in Fig. [2j The crossover between the single- and multi-peak states 
in the plane of (a\, /), which summarizes results of many simulations performed at different 
values of the parameters, is displayed in Fig. [3j 

In fact, the transition to the multi-peak patterns with the increase of the amplitude is 
a straightforward manifestation of the modulational instability of extended states, which is 
a property well known both in conservative models and in those based on the CGL 
equations [3], provided that the nonlinear terms in the model have the sign corresponding 
to the self-focusing. As mentioned above, this is indeed the case for the self-interaction of 
magnons. 
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FIG. 2: (Color online) Normalized profiles of stable localized states, obtained with four different 
values of the strength of the pump, /, in the model with a broad source, corresponding to A = 80 
in Eq. (|4|). Other parameters are the same as in Fig. Q3 Also shown, for comparison, is the shape 
of the source S(z). The transition from the single-peak profile to a multi-peak one occurs at some 
point in interval 0.05 < f < 0.10. 




FIG. 3: Stationary localized states generated by the simulations of Eqs. ([2]) and ([3]) feature single- 
and multi-peak shapes in areas I and II, respectively. Other parameters are fiQ = 0, rj = 0.007, 
r = 0.22, and a 2 = 3. 

B. Asymmetric solutions 

The above analysis of the phenomenological model predicted a possibility of the tran- 
sition from symmetric to asymmetric states, at least in the strongly dissipative system. 
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FIG. 4: The border between the areas of monost ability (I) and bistability (II), in the plane of 
parameters (o"2,o"i)- Other coefficients are = 0, 77 = 0.007, r = 0.22, / = 0.15, and A = 80. 
In area II, asymmetric initial conditions, with one component equal to zero, generate persistently 
asymmetric nonstationary localized solutions, while in area I the same initial conditions converge 
to stationary symmetric states. In either area, symmetric initial conditions always converge to a 
stationary symmetric localized state, which may feature a single- or multi-peak shape in area I, 
and is always of the multi-peak type in area II. 



This possibility was systematically tested by running simulations of Eqs. ([2]) and (J3j) with 
asymmetric initial conditions [in fact, this was done by choosing = 0) = 0]. We 

have found that, in some region of the parameter space, any initial condition, symmetric or 
asymmetric, converges to a single symmetric localized state. On the other hand, another 
parameter region exists too, where the asymmetric initial configuration generates solutions 
which remain asymmetric indefinitely long. In fact, this situation implies a bistability in 
the system, as the symmetric input would always generate a stationary symmetric state, 
at the same values of the parameters. The phase diagram in the parameter space (<j 2 ,cri) 
showing the regions of the monostability and bistability is displayed in Fig. HI It is relevant 
to stress that, with the variation of the parameters, the transition from single- to multi-peak 
symmetric states in the class of symmetric solutions (see above) always precedes the onset of 
the bistability. In other words, in the case of the bistability the symmetric states coexisting 
with asymmetric ones always represent multi-peak patterns. 

In fact, the persistently asymmetric states, unlike the symmetric ones, never relax to a 
stationary shape. While keeping an overall localized form, they demonstrate quasi-regular 
oscillations between multi-peaked configurations in the two components, and Typ- 
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FIG. 5: (Color online) Typical examples of the dynamical behavior of asymmetric states, generated 
by initial conditions with ^-(z, t = 0) = . Parts (a,b) and (c,d) pertain to the pump's strength 
/ = 0.09 and 0.12, respectively. Other parameters are fio = 0, r\ = 0.007, r = 0.22, a% = a? = 3.5, 
A = 80 . Panels (a) and (c) display the spatiotemporal evolution of absolute values of the two 
fields, | ^±(2, t)\, while panels (d) and (d) show the profiles of the solutions observed at t = 50, 000. 

ical examples of the dynamical behavior of asymmetric states are displayed in Fig. [5j It is 
worthy to note that, at smaller values of the pumping strength, /, each component *f>± of 
the asymmetric state generates traveling waves, see Fig. [5](a), while at larger values of / the 
nonstationary dynamics is represented by standing waves, as seen in Fig. EJc). 

In addition to the asymmetry between the components and it is relevant to 
monitor the spatial symmetry of these states. The conclusion is that the solutions are 
spatially asymmetric at lower values of the forcing, 0.07 < f < 0.18, and they become 
spatially symmetric, with \P±(— 2) = ty±(z), at higher values of /, although the solution 
remains nonstationary, with unequal components and Pairs of panels (a,c) and 
(b,d) in Fig. [5] represent, as a matter of fact, typical examples of the spatial asymmetry 
and symmetry, respectively (in other words, traveling-wave and standing-wave patterns, as 
mentioned above). 
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V. COMPARISON WITH THE EXPERIMENT 



The model presented above can be used for the description of recent experiments demon- 
strating spatiotemporal patterns in magnon condensates formed through the parametric 
pumping Jjj]. The respective pumping field was created by microstrip resonators of two 
types with different spatial profiles of the microwave field. The spatiotemporal evolution 
of the condensates was studied, using a pumping pulse of a finite width and applying the 
technique of the time- and space-resolved Brillouin light scattering. 

Before comparing the theoretical results with experimental data, it is necessary to stress 
that the source term in Eqs. (J2J) and (jl]), fS(z), does not directly represent the pumping 
field of the resonator used in the experiment. There are two reasons for that. First, in the 
experiment the resonator excites pairs of primary magnons by means of the parametric 
pumping. The process of parametric pumping has a threshold, with respect to the pumping 
field. Thus, although, as discussed above, the condensate can be created at any value of / 
which determines the strength of the source in Eqs. ([2j) and ([3]) , from the experimental point 
of view the process has a threshold with respect to the microwave pumping field. Second, 
the primary magnons create the condensate as a result of the multi-step thermalization 
process [? ], which is accompanied by the spreading of the magnon cloud in the physical 
space. Strictly speaking, to calculate the spatial profile of the pumping, one needs, first, 
to solve the system of equations similar to Eqs. (|2J) and which describe the excitation 
of the primary (parametric) magnons. Then, the so generated profiles should be used in 
the model accounting for the thermalization process. In its entire form, this problem can 
hardly be solved even numerically. However, since the group velocity of magnons can be 
easily found from their spectrum, and the thermalization time is measured experimentally 



29j , the spreading of the magnon cloud during the thermalization can be estimated as being 
below 100 /im. 

The presence of the effective threshold, can be easily taken into account. In Ref. [h]] it 
was shown that the number of magnons in the condensate depends on the strength of the 



microwave field, h, as \Jh 2 — /z 2 hr . Accordingly, we can represent the source term as: 



fS{z) ~ y/h'{z) - (O 2 , (26) 
where the effect of the spreading of the magnon cloud during the thermalization process is 
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taken into account by the introduction of the effective threshold field . For values of 
z at which h 2 (z) — (h^. ) 2 is negative, Eq. fl26|) is replaced by S(z) = 0, since no primary 
magnons are excited in this region. 

As mentioned above, the parametric pumping excites pairs of primary magnons with 
wavenumbers +ko and — ko, which makes experimentally relevant initial conditions for the 
creation of the condensate symmetric. Therefore, it is not surprising that only symmetric 
solutions have been observed in the experiment. 

For the sake of the qualitative comparison, one needs to determine parameters in Eqs. ([2]) 
and ([3]). The coefficients of the linear and nonlinear dissipation, rj and r, can be extracted 
from the experimental data via the analysis of the decay of the condensate when the source 
(pumping field) is switched off 22]. For the YIG films used in the experiments, these 
coefficients are found to be 77 ~ 0.007 and r ~ 0.22 (in the scaled notation, corresponding 
to the characteristic nonlinear time of r = 3.6 ns). Finally, as mentioned above, due to 
the fact that the nonlinearity in the system is caused by the four-magnon interaction (cubic 
nonlinearity) , the XPM/SPM ratios are o\ = 02 = 2. 

The comparison between the predictions of the phenomenological theory and the experi- 
mental data borrowed from Ref. [3] is illustrated in Fig. [6J which shows the spatial profile 



of the stationary symmetric state calculated using the above parameters and expression ([26]) 
for the source. The distribution of the parallel component of the microwave field of the wire 
resonator was used as h(z), see Fig. 3(a) in Ref. [12|. The experimentally measured profile 
of the condensate density is also shown in in Fig. [6l Threshold in Eq. (|26l) was treated 
fitting parameter. 

Since Eqs. (|2]) an d (G2) describe the spatiotemporal dynamics of the condensate, which 
has been also been explored experimentally in Ref. [l2], it makes sense to compare spa- 
tiotemporal patterns too. Figure [7J represents examples of the spatiotemporal evolution of 
the magnon field, as predicted by Eqs. ([2]) and ([3]), and the respective experimental data. 
The situation presented in the figure corresponds to the case of the excitation of magnons 
by the microstrip resonator [l^]. Accordin gly, the distribution of the microwave field cor- 
responding to that shown in Fig. 3(b) of 12j was used for the calculation of the source 
term fl26|) . Point t = in Fig. [7J corresponds to the start of the pumping pulse. The 
delay in the emergence of the condensate is caused by the process of the thermalization of 
magnons, through their nonlinear interactions, as discussed above [2]. After about 1000 ns, 
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FIG. 6: A typical example of the comparison of the normalized stationary profile of the magnon 
field, as predicted by the phenomenological equations, and its experimentally observed counterpart. 
Further details are given in the text. 

one observes decay of the condensate density, which is due to the fact that the pumping is 
switched off at that time. While calculating the theoretical profiles, we have concluded that 
the most appropriate results are produced by the simulations with 5fiQ = —0.011, which 
corresponds to —3 x 10 6 s _1 , in the physical units. This value is in agreement with that fact 
that the pulsed character of the experiment, with characteristic time t = 1 /is, results in the 
uncertainty of magnon states within a frequency interval Af = 2n /t ~ 6.28 x 10 6 s -1 . 

Comparing the theoretical results and experimental data displayed in Fig. [3, one can 
conclude that the present model describes the spatiotemporal evolution of magnon conden- 
sate fairly well. It is clearly seen that the spatial size of the condensate increases with the 
growth of the pumping power. This property is explained by the fact that the spatially 
uniform enhancement of the pumping field, just by increasing the microwave power sent to 
the resonator, changes the actual profile of the corresponding source. Indeed, for a larger 
pumping power, condition h(z) > see Eq. f )26|) . is satisfied in a larger interval of z. 

VI. CONCLUSION 

In this work we have proposed a system of semi-phenomenological one- dimensional equa- 
tions for the dynamical description of the formation of spatiotemporal patterns in Bose- 
Einstein condensates of magnons. The equations combine essential features of the GPE- 
and CGL-type models. The equations include localized source terms, which represent the 
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FIG. 7: (Color online) Spatio-temporal evolution of the density of the magnon condensate in 
the case of the pulsed excitation of the microstrip resonator by the microwave field, for different 
pumping powers P p , as indicated in the figure. The left and right columns display the experimental 
and theoretical results, respectively. Further explanations are given in the text. 

microwave field pumping magnons into the condensate. In the limit case of the source 
represented by the ^-function, we have found exact analytical solutions for symmetric lo- 
calized states of the condensate. The modes also predicts the possibility of the existence 
of asymmetric solutions, with unequal amplitudes of the left- and right-traveling magnon 
waves. Systematic simulations of the model have demonstrated that the analytically found 
solutions are always stable. 

Replacing the 5-function source by the broader one, which is relevant to the experimental 
situation, we have found a transition from the single-peak stationary symmetric solution to 
stationary multi-peak patterns, which may be explained as a manifestation of the mod- 
ulational instability of broad nonlinear modes. While, in direct simultations, symmetric 
initial conditions always converge into symmetric stationary states, in a part of the param- 
eter space asymmetric inputs may generate persistently asymmetric nonstationary localized 
modes, which may be realized as patches of traveling or standing waves, in the cases of low 
and high amplitudes, respectively. 

Comparison with experimental observations demonstrates that the spatial and spatiotem- 
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poral patterns predicted by the phenomenological equations may provide for a good model 
of the experiment, with some parameters found from fundamental characteristics of the fer- 
romagnetic medium, and some others used for the fitting. Moreover, the model suggests to 
look for new types of patterns in the experiment, such as asymmetric ones. 
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